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Abstract: A holographic dual description of a 2+1 dimensional system of strongly 
interacting fermions at low temperature and finite charge density is given in terms of 
an electron cloud suspended over the horizon of a charged black hole in asymptotically 
AdS spacetime. The electron star of Hartnoll and Tavanfar is recovered in the limit 
of zero temperature, while at higher temperatures the fraction of charge carried by 
the electron cloud is reduced and at a critical temperature there is a third order 
phase transition to a configuration with only a charged black hole. The geometric 
structure implies that finite temperature transport coefficients, including the AC 
electrical conductivity, only receive contributions from bulk fermions within a finite 
band in the radial direction. 
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1. Introduction 

There has been considerable recent interest in developing holographic models of 
strongly coupled physics in low dimensions with a view towards condensed mat- 
ter systems (see [1, 2, 3, 4] for reviews). This can be motivated both from the point 
of view of extending the gauge theory /gravity correspondence to include a variety of 
interesting field theoretic systems without supersymmetry, and also from the point of 
view of gaining a new theoretical handle on materials containing strongly correlated 
electrons. 

In a recent paper Hartnoll and Tavanfar [5] considered a simple holographic 
model for strongly interacting fermions in 2+1 dimensions at zero temperature and 
finite charge density. In their approach, which builds on earlier work in [6, 7], the 
bulk Maxwell field, which is dual to the field theory current, is sourced by an ideal 
fluid consisting of charged free fermions. The combined Einstein, Maxwell, and fluid 
field equations in the bulk spacetime have planar solutions, referred to as electron 
stars in [5] , where the charge and energy densities of the bulk fermion fluid have a 
non-trivial radial profile. The geometry is asymptotically AdS but deep in the elec- 
tron star interior the metric exhibits Lifshitz scaling with a non-universal dynamical 
critical exponent that depends on the couplings of the model. Similar constructions 
were considered in [8, 9] for neutral and charged free fermion fluids supported by a 
degeneracy pressure. 

In the present paper we extend the gravity dual description of [5] to include 
finite temperature configurations in the boundary field theory. We construct static 
solutions of the bulk field equations where an 'electron cloud' is suspended above 
the horizon of a charged black hole, or more precisely a black brane with a planar 
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horizon. The electron cloud has both an outer and an inner edge. The outer edge is 
also found in the electron stars of [5] but the inner edge is a new feature, found only 
at finite temperature. At the inner edge the gravitational pull of the black brane on 
the electron fluid is balanced by electrostatic repulsion. 

In the fluid description, observables in the boundary theory, such as the electric 
conductivity, only receive contributions from bulk fermions located within a band 
of finite width in the radial direction. The sharpness of the edges is presumably 
an artifact of the classical perfect fluid description and we expect both quantum 
corrections and fluid interactions to give rise to tails in the bulk fermion profile that 
fall off towards the boundary and the black brane horizon, respectively. 1 

At the level of classical geometry, a finite temperature in the boundary theory 
is introduced by including a black hole in the bulk spacetime, with a non-vanishing 
Hawking temperature. Quantum effects in the bulk include thermal Hawking radi- 
ation, which comes to equilibrium with the bulk fermion fluid, but at weak gravi- 
tational coupling this thermalization in the bulk is suppressed and will not be con- 
sidered here. This simplifies our analysis considerably as it allows us to use a zero 
temperature equation of state for the free fermions in the bulk to capture finite 
temperature effects in the boundary field theory. 

At low temperatures in the boundary theory, the electric charge in the bulk 
geometry is partly carried by the electron cloud and partly by the charged black 
brane inside it. As the temperature is raised, the two edges of the electron cloud 
move towards each other and an ever larger fraction of the total electric charge in 
the bulk geometry resides inside the black brane. The two edges of the electron 
fluid meet at a finite critical temperature, above which there is only a black brane 
solution with no electron cloud present. We find that the system undergoes a third 
order phase transition at the critical point. 2 

As the temperature is lowered, on the other hand, the inner edge of the electron 
cloud approaches the black brane horizon, which in turn recedes towards vanishing 
area. In the zero temperature limit, one recovers the electron star geometry where 
there is no longer any black hole and all the charge is carried by the electron fluid. We 
confirm this by showing that the horizon recedes from the boundary, the geometry at 
low temperature encodes the dynamical exponent z, and by comparing free energy 
densities. We find that at low temperatures the free energy density of an electron 
cloud geometry smoothly goes over to that of an electron star. Furthermore, we 
compare the free energy densities between an electron cloud solution and an AdS- 
RN black brane solution at finite temperature. We find that, whenever a solution 
with an electron cloud exists, it is favored over an AdS-RN black brane. 

1 Such tails are, for instance, found in an alternative approach to including bulk fermions based 
on a single fermion wave equation [11]. 

2 In an earlier preprint of this paper it was incorrectly stated that the phase transition was second 
order. The correct behavior was identified in [10]. 
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The electrical conductivity at finite temperature can be obtained for this system 
using, by now, standard holographic techniques (see for instance [1, 2]). We find that 
the finite temperature conductivity smoothly interpolates between the AdS-RN and 
electron star results [5]. 

The same finite temperature solutions were found independently by Hartnoll and 
Petrov in [10]. Initially there was a discrepancy between their work and ours in the 
analysis of the phase transition, but after correcting an error in our expression for 
the free energy density, we now also find a third order phase transition. 



2. Field equations and electron cloud solutions 

The Einstein-Maxwell equations with a negative cosmological constant and a charged 
perfect fluid are 

R,u - \g, v R - ^ = K\Tfr^ X + T*: id ) , Vi^ = e 2 Jf id . (2.1) 

We adopt units where the characteristic AdS length scale is L = 1. The source terms 
are given by 

^Maxwell = i_ (FMAF / _ lg^ Fx ^) , (2.2) 

T^ d = (p + p)u li u v +pg IM/ , (2.3) 
J** = au, , (2.4) 

where a is the charge density of the fluid, p is its energy density, p the pressure, and 
the four velocity, u^u, — — 1. The justification and limitations of the perfect fluid 
description are discussed in detail in [5] and the same considerations apply here. 
We look for static black brane solutions with planar symmetry, 

ds 2 = -f(v)dt 2 + g(v)dv 2 + \{dx 2 + dy 2 ) , A = -h(v) dt , (2.5) 

where the radial coordinate goes from v — > at the asymptotic boundary to a 
constant value v — vq at the black brane horizon. We find it convenient to introduce 
a scale invariant variable u — — log(i>/i>o), such that u — at the horizon and u — > oo 
at the boundary, and work with rescaled fields, 

/ = v%f, g = v 2 g, h = v h, p = K 2 p, p = K 2 p, a = ena . (2.6) 

The equations of motion (2.1) can then be expressed in a first order form, convenient 
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for numerical evaluation, 



' • ^ • /u ■■>» 2 "<n < 2 \fm>. (2.7) 



du 2 



dk , r. -2u f l lt , A 9& 



du + k ' 2 "U<><- • ( 2 - 8 ) 



lf + l -f-A = e-^ ) (2.9) 
y du g du I j 

where k = dh/du. Following [5], we assume a free fermion equation of state defined 
via 

a = (3 deey/e 2 — rh 2 , P = P dee 2 \ / e 2 — rh 2 , —p = p — jia, (2-10) 



where /3 is a coupling dependent dimensionless constant, m is proportional to the 
electron mass, m 2 = ^m 2 , and the (rescaled) local chemical potential fi is given by 

the background Maxwell gauge field in the tangent frame, fi = h/ \J~f . 
As discussed in [5], there is a range of parameters, 

e 2 ~y<l, (2.11) 

for which we can assume a classical bulk geometry with a non-trivial back-reaction 
due to the fermion fluid. If, at the same time, the Compton wavelength of the 
fermions is small compared to the AdS length scale, 

mi»l, (2.12) 

then we are also justified in taking spacetime to be locally flat in the fermion equation 
of state. These conditions amount to the dimensionless parameters in the equation 
of state (2.10) taking order one values [5], 

/3~1, m 2 ~l. (2.13) 

The construction of the electron cloud geometry proceeds in a few steps. First 
we solve the vacuum equations, with a = p = p = 0, to find the charged AdS-RN 
black brane solution inside the cloud, 



/ = e 2M +^e- 2u -(l + ^)e-", g=C, h = q(l - e' u ) . (2.14) 



The dimensionless constant q is proportional to the charge carried by the black 
brane and we have used the freedom to rescale the time coordinate t to fix the 
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Figure 1: The auxiliary function r(u) in (2.16), used in the construction of the electron 
cloud solution, is plotted for rh = 1 (dashed red) and rh = 0.55 (solid blue). 



overall normalization of /. When the charge parameter is in the range q 2 < 6 the 
black brane is non-extremal with a non-degenerate horizon, where the local chemical 
potential p, vanishes. As we move away from the horizon the chemical potential grows 
but remains too small to support a fermion fluid until 

/i 2 = S^>m 2 . (2.15) 

m 

We will only consider non-vanishing fermion mass. For zero fermion mass the inner 
edge of the electron cloud reaches the horizon for all temperatures. The geometry, 
with the back-reaction from the fermion fluid included, is then harder to determine 
and we will not consider this case here. 

The condition (2.15) is easily seen to be equivalent to 

„ 2 m 2 e u (e 2u + e u + 1) 

q > - 9 - = r(u). (2.16) 

The right hand side is shown for two different values of rh 2 in Figure 1. It can be 
read off from the figure that: 

• There cannot be any fermion fluid outside a non-extremal black brane if rh 2 > 1, 
since in this case r(u) > 6 for all u > 0. This restriction on rh was already 
seen in [5] as a condition for the existence of electron star solutions at zero 
temperature. 

• For rh 2 < 1 and a near-extremal black brane with q 2 < 6, the condition (2.16) 
is satisfied within a finite interval u\ < u < U2, indicated in the figure. The 
endpoints of the interval correspond to the inner and outer edges of the electron 
cloud in a "probe" approximation, where the back-reaction on the geometry 
due to the fermion fluid is ignored. 
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Figure 2: The radial profiles of the fluid variables (a,p,p) for rh 
q 2 = 4.49. 



0.55, j3 = 10, and 



• For rh 2 < 1 and q 2 below a critical value (which depends on m 2 ), the condition 
(2.16) is not satisfied for any u > 0. This implies there is a critical temperature 
above which there is only a black brane and no electron fluid in the bulk. 

The second step in the construction of an electron cloud solution with back- 
reaction included is to numerically integrate the field equations (2.7) - (2.9) starting 
from the inner edge of the electron cloud. The functions /, g, h, and k are continuous 
at the matching point u = U\ and thus we can generate initial values for the numerical 
integration using the exact AdS-RN solution (2.14) with q determined from (2.16) 
evaluated at u = U\. 

The local chemical potential fx goes to zero in the asymptotic u — > oo region and 
the numerical integration is terminated at a point u = u s where the condition (2.16) 
is no longer satisfied. We find that the back-reaction of the fermion fluid on the 
geometry leads to u s > u 2 , and that this effect becomes more pronounced, u s ^> u 2 , 
at low temperature. Figure 2 shows numerical results for the fluid variables a, p, 
and p for rh = 0.55, (3 = 10, and q 2 = 4.49. 

The third and final step in the construction is to obtain the spacetime geometry 
outside the electron cloud by matching the numerical solution onto a charged black 
brane solution at u = u s in much the same way as is done for electron stars in [5]. 
The exterior solution has the general form 



./ 



2u 



9 



c 2 e 4u 



f 



h = /x s - Q s e~ 



(2.17) 



where the subscript s on the various constant parameters is a reminder that they 
are determined by matching onto a numerical solution at u = u s . We already fixed 
the overall scale of the time coordinate when writing the inside black brane solution 
in (2.14) so now there appears an extra parameter c s in /. Also, since the external 
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solution only extends to u — u s and not to an event horizon, we do not require the 
usual relationship between M s , Q s , and /a s found for a vacuum AdS-RN black brane. 
The parameters in (2.17) are instead determined to be 

<? a = f(u 8 )9(u a )e-^ , (2.18) 
Q s = k{u s )e u % (2.19) 
Us = h(u s ) + k(u s ) , (2.20) 

M s = f(u s )g(u s )e- u ° + \h\u a )e u - - f(u s )e u ° . (2.21) 

These parameters refer to the rescaled fields in (2.6) while the physical parame- 
ters appearing in an external AdS-RN solution with a canonically normalized time 
coordinate are given by 

H= J ±, Q = — 2 , M=^- 3 . (2.22) 
c s v c s v£ cjv$ 

Once the parameters of the external black brane solution have been determined for 
given values of m, (3, and q, standard methods can be used to obtain the free energy 
density as a function of temperature for these geometries. This will be carried out 
in Section 3 below. 

The next step is to determine the Hawking temperature of the electron cloud 
geometry, which is to be identified with the temperature in the boundary field theory. 
The Hawking temperature is easily obtained from the behavior of the Euclidean 
metric near the horizon. One finds 

T -^, (2.23) 

where we have again divided by the physical chemical potential /i in order to have a 
dimensionless quantity to work with. 

In the limit of zero temperature we expect to recover the electron star geometry 
of [5]. In this case the radial coordinate extends to u — > — oo and the metric exhibits 
Lifshitz scaling in the deep interior. 

In Figure 3 one can see that the curvature scalar R as a function of the proper 
distance s measured from the outside of the electron cloud is approaching the ex- 
pected asymptotic value for a Lifshitz geometry as we lower the temperature T before 
plunging to R = —12 at the horizon. From R evaluated on the solution and the met- 
ric, one confirms the correct dynamical exponent, e.g. z = 5.75466 when rh = .55 
and $ = 10. 

On the left in Figure 4, we see that the horizon radius vanishes as the temperature 
is lowered indicating that the horizon recedes from the AdS boundary. Furthermore 
by fitting the low temperature data in the figure, we find that 

-ocf— V, (2.24) 
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Figure 3: The curvature scalar R versus the proper distance s measured from the outer 
edge of the electron cloud for rh = .55 ,f3 = 10 and for T/fi values .55, .22, 9.4 x 1CT 2 , 9x 1CT 4 
and 1 x 10" 5 . The curves extend further to the right with decreasing temperature. The 
value for R in the Lifshitz region deep inside an electron star with the same rh and /3 is 
shown as a horizontal line for reference. 
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Figure 4: On the left we plot the horizon radius Vq as a function of the temperature, T. 
On the right we show the charge parameter of the inside black brane solution (2.14), q 2 , 
versus T for rh = .55 and (3 = 10. For reference we plot the extremal value q 2 = 6. 



where z matches the appropriate expected Lifshitz exponent to a very high precision. 
This is further evidence that the anisotropic scaling found inside the electron star is 
recovered by our electron cloud solutions at low temperature. 

The plot on the right in Figure 4 shows that the zero temperature limit in 
the boundary theory is in fact not obtained by approaching an extremal interior 
black brane, which would have q 2 = 6. This is at first sight counterintuitive but can 
ultimately be traced to the focusing effect that the electron fluid has on the geometry 
outside the black brane. If one were to try to construct an electron cloud solution 
starting with a value of q 2 that is closer to 6, than what is seen in the figure the 
metric would collapse to a curvature singularity at a finite proper distance outside 
the horizon and the solution would never reach an asymptotic AdS region. 
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3. Free energy 



Further evidence that the electron cloud solution is the proper finite temperature 
extension of an electron star comes from comparing free energy densities. We obtain 
the free energy by evaluating the on-shell Euclidean action of the bulk system, in- 
cluding the usual Gibbons-Hawking boundary term [12] and boundary counterterms 
required for regularization [13, 14]. A bulk action for the charged electron fluid also 
needs to be included, as described in [5]. The full bulk action turns out to be the 
integral of a total derivative, and when combined with the appropriate boundary 
terms, it gives a simple result for the free energy density, 

F = M -fxQ -sT, (3.1) 

where s is the Bekenstein-Hawking entropy density. 3 This can be simplified by using 
the thermodynamic relation 



3 M - fiQ - sT = 0, (3.2) 



2 

giving 

M , . 

F = -~. (3.3) 

The relation (3.1) follows from the radial conservation of the quantity 

D = -^L (-2hk - 2/ + . (3.4) 

By using the equations of motion (2.7) - (2.10), it is straightforward to check that 
= and one then evaluates D at the horizon and at the u — > oo boundary to 
obtain (3.2). 

Using (2.22), the free energy density can be re-expressed in terms of output 
parameters from our numerical evaluation, 



1 M s 

F = ~^^A- 3.5 

2 C y 



The factor of Vq in the denominator tells us that we should instead work with the 
dimensionless quantity 

1 c s M s 

(3.6) 



IM* 2 ^ > 

when comparing free energy densities. 

In Figure 5 we use these dimensionless variables to compare the free energy 
densities of various geometries for a typical case when m = .55 and (3 = 10, holding 

3 In an earlier version of the paper the sT term was missing from the expression for the free 
energy density. This led to an incorrect characterization of the phase transition. 
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Figure 5: The free energy densities of AdS-RN black brane and electron cloud solutions 
for = 10 and rh = .55. In addition the free energy for an extremal black hole and the 
electron star solution of [5] are shown with a red box and a green dot, respectively. 



the chemical potential fixed. One readily sees that the electron cloud solution is 
preferred over the black brane solution up to the point where the local chemical 
potential is too low to support any fluid. Beyond this point the only solution is an 
AdS-RN black brane. At low temperatures, on the other hand, the free energy density 
of the electron cloud geometries approaches that of the corresponding electron star. 

In addition to the low-temperature regime, it is also interesting to ask about the 
nature of the transition to the AdS-RN black brane solution at higher temperatures. 
To address this issue, we consider the difference in energy densities between an 
electron cloud solution just below the critical temperature T c and an AdS-RN black 
brane solution, 

(3.7) 



/i 3 / V^ 3 / AdS-RN V^ 3 / EC 

at the same value of T/ri. Figure 6 shows a log-log plot of this difference near the 
critical point where one loses the cloud solution at T c //x = 0.058868 for rh = .55 and 
(3 = 10. The solid curve in Figure 6 is a straight line of slope 3 giving numerical 
evidence of a third order phase transition where 

1^— H\ (3.8) 



/i 3 y v (J- 



This feature was observed numerically in [10], and those authors also gave a simple 
analytic argument for this behavior. 
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Figure 6: The difference between the black brane and electron cloud free energies near 
the phase transition temperature at (3 = 10 and rh = .55. 



4. Electric conductivity 



The finite temperature AC conductivity at zero momentum can be computed in an 
analogous way as was done for the zero temperature electron star in [5]. In fact, 
the finite temperature computation is more standard since in this case the ingoing 
boundary conditions [15] for the fluctuations in the gauge field are imposed at a 
smooth black brane horizon rather than at the (mildly) singular Lifshitz horizon 
inside an electron star. 

The background is perturbed, assuming a time dependence of the form e~ luJt , 
and the resulting equations are linearized. To get a closed system of equations, the 
following perturbations are needed 



A x = 8AJu)e 



-iwt 



9tx = 5g tx (u)e 



-iwt 



u x = Su x (u)e 



-iuit 



This leads to a system of four first order differential equations 

5A X + j5g tx + e 2u fx5u x 

dSgtx Q ~ lo^jcJ 

— 25g tx + 2 — 6A X 

du du 




d5B x 
du 



+ 



2e 



2u 



dh 
du 



-u 2 s g 



+ 



h 



5A X = 0, 



(4.1) 

(4.2) 

(4.3) 
(4.4) 

(4.5) 



where u s is defined in terms of the canonical normalized frequency u as u s = c s VqU. 
Equation (4.4) can be regarded as the definition of the auxiliary function 5B X . We 
note that (4.4) and (4.5) form a closed system involving only 5A X and 5B X . 
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Figure 7: Real and Imaginary part of the conductivity for rh = 0.55, (3 = 10. The curves 
blue dots, dashed green and solid red correspond to temperature value T/Tc = 2, 1 and 
0.162, respectively. Tq/h = 0.05887 denotes the critical temperature where the phase 
transition occurs. Curves for lower values of T are almost indistinguishable from the solid 
red curve. 



At the horizon, ingoing boundary conditions are imposed [15]. This implies 
5A X — » u^^r and 5B X — > iuvr^r as u — > and T being the Hawking tempera- 
ture of the AdS-RN black brane solution (2.14). At the AdS boundary, where the 
background is of the form (2.17), the behavior of those functions is 

5A X = if + A x ^e~ u + • • • , (4.6) 
SB X = fl(°> + B^e~ u + ■■■ . (4.7) 

The coefficients A% and B x are connected, e.g. Bx = c s A x l \ This relation can be 
used to express the conductivity as 

1 JDX I A o\ 

which is manifestly invariant under the rescaling described in the previous sections. 

Our results for the conductivity are obtained by numerics. This is achieved by 
integrating out from the horizon in the background of an AdS-RN solution, as already 
indicated, until the inner edge of the electron shell is reached. There, 5A X and 5B X 
need to be continued smoothly into a solution of (4.4) and (4.5) with the electron 
cloud solution as background. At the outer edge, a second matching to the exterior 
solution must occur. Finally, the coefficients A x ^ and B x ^ can be read off at the 
boundary and plugged in to (4.8). 

A plot of the conductivity can be seen in Figure 7. The pole in the imaginary 
part, as usual, indicates the presence of a delta peak in the real part. The offset 
in the conductivity goes rather quickly to zero once the electron cloud is in place. 
This is also shown in Figure 8. Parameterizing the real and imaginary part of the 



- 12 - 




1(T 6 10- 5 ]0~ 4 0.001 0.01 




Figure 8: Numerical results for the coefficients <to, °"2 an d ct-% in the electrical conductivity 
in equation (4.9). 

conductivity as 



for small u, it can be seen from our numerical data that 02 and <r_i level off for 
small temperatures and o"o decreases like T 2 in this limit. This is consistent with the 
zero-temperature result in [5], where <To is not present, and with u/T scaling at low 
temperatures. Above the critical temperature the background is an AdS-RN black 
brane and the calculation of the conductivity reduces to the one described in [2]. 

5. Discussion 

One of the most interesting recent developments in holographic model building is the 
observation that at low temperatures it can be energetically favorable for black holes 
in AdS space to eject their charge in the form of matter "hair" [7, 16, 17, 18, 19]. In 
the boundary theory this hair may give rise to many interesting features including 
spontaneous symmetry breaking [20] , dynamical critical exponents [5] and non-fermi 
liquids [7, 21]. In general, this may allow one to use holographic techniques to study 
condensed matter systems not amenable to other theoretical tools. 

In this paper we explored the finite boundary temperature generalization of the 
electron star configuration described in [5]. The electron star solution is a zero 
temperature model for a quantum phase transition displaying dynamical critical 
exponents as well as non-fermi liquid features. In the bulk the configuration is that 
of a zero temperature ideal charged fermion fluid. In the deep interior this fluid has 
an asymptotic Lifshitz scaling symmetry. As one approaches the boundary, the gauge 
potential is screened by the charged fluid. Eventually, the local chemical potential 
falls below the fermion mass causing the fluid density to vanish. In the asymptotic 
region one is left with an AdS-RN geometry. 

At finite boundary temperatures there is instead a cloud-like configuration, with 
the electron fluid hovering outside an AdS-RN black brane geometry. This configu- 
ration is only possible when the fluid and black brane have same sign charges such 





(4.9) 
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that electrostatic repulsion balances gravitational attraction. Within the fluid and 
its exterior, the electron cloud solution is similar to the electron star. The gauge 
field is screened and eventually the fluid can no longer be supported. We found that 
as one lowers the temperature, for fixed chemical potential one smoothly obtains the 
free energy of the electron star solution at zero temperature. 

In the other extreme, beyond a critical temperature the local chemical potential 
is always lower than the fermion mass and no fermion fluid is supported. In this 
case we are left with an AdS-RN black brane geometry. We studied the transition 
numerically and found that it is a third order phase transition, as pointed out in [10]. 

We can summarize the phase diagram as follows, in the high temperature regime 
there is an AdS-RN black brane. As one lowers the temperature it is favorable for 
the black brane to expel some of its charge in the form of an electron cloud hovering 
over the horizon. As one cools the system further, the interior black brane carries 
less and less charge and has shrinking area. Finally, at zero temperature the black 
brane is gone and the fluid takes the form of an electron star. 

In addition to studying electron cloud thermodynamics we also computed the 
conductivity, finding that it is nicely consistent with the electron star results of [5] 
at low temperature and with the AdS-RN black brane at high temperature. 

In [21] it was argued that if the fluid in an electron star experiences a local 
temperature and magnetic field, it is possible to detect a Fermi surface evidenced by 
Kosevich-Lifshitz oscillations. It should also be possible to see such oscillations in 
electron clouds. 

A major challenge for the construction we are working with is how to interpret 
the electron fluid directly in terms of operators in the boundary field theory. It would 
be interesting to see how the work of [8] may be extended for charged fluids in order 
to further characterize the nature of the underlying quantum critical point. 
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